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Abstract 

The chaotic behavior of the modified Rayleigh-Duffing oscillator with 
</> 6 potential and external excitation which modeles ship rolling motions are 
investigated both analytically and numerically. Melnikov method is applied 
and the conditions for the existence of homoclinic and heteroclinic chaos are 
obtained. The effects of nonlinear damping on roll motion of ships are an¬ 
alyzed in detail. As it is known, nonlinear roll damping is a very important 
parameter in estimating ship reponses. The predictions are tested numeri¬ 
cal simulations based on the basin of attraction. We conclude that certains 
quadratic damping effects are contrary to cubic damping effect. 

keywords Modified Rayleigh-Duffing oscillator, Melnikov criterion, nonlin¬ 
ear damping effect, ship rolling, basin of attraction 

1 Introduction 

The analysis of nonlinear dynamic systems is now a major theme in both aca¬ 
demic and industrial perspective and touches many areas, such as hydrodynamics, 
aerospace, civil engineering, transport, musical acoustics, nuclear engineering and 
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others [Enjieu Kadji and Nana Nbendjo, 2012; Hayashi, 1964; Nana Nbendjo et 
al., 2007; Nayfeh and Mook, 1979; Tchoukuegno, 2002; Tchoukuegno, 2003; 
Yamapi, 2003; Yamapi, 2007], The first works date from the nineteenth century 
including Poincare, but currently it is experiencing a resurgence of interest due to 
the need to optimize, streamline structures commonly used and subjected to sig¬ 
nificant levels of excitement, or control the instabilities oscillations. Modeling and 
study of the behavior of a physical system is a major problem even complicated in 
case the present system nonlinearities. Much of the discussion in the physics and 
engineering literature concerning damped oscillations, focuses on systems subject 
to viscous damping, that is, damping proportional to the velocity, even though 
viscous damping occurs rarely in physical systems. Phenomenological models 
describing some type of nonlinear dissipation have been used in some applied sci¬ 
ences such as ship dynamics [Bikdash et al., 1994; Falzarano et al., 1992], where a 
particular interest has deserved the role played by different damping mechanisms 
in the formulation of ship stability criteria, and vibration engineering [Ravindra 
and Mallik, 1994a; Ravindra and Mallik, 1994b]. Damping in certain applied sys¬ 
tems plays an important role, since it may be used to suppress large amplitude 
oscillations or various instabilities, and it can be also used 

as a control mechanism [Litak et al., 2009; Miwadinou et al., 2015; Rand et 
al., 2000; Sanjuan, 1999; Soliman and Thompson, 1992, Taylan, 2000]. 

This paper is concerned with the appearance of homoclinic and heteroclinic 
instabilities and chaos in a triple-well modified Rayleigh-Duffing oscillator. Many 
problems in physics, chemistry, biology, Engineering, etc., are related Rayleigh- 
Duffing oscillator. [Miwadinou et al., 2014; Miwadinou et al., 2015; Siewe Siewe 
et al., 2009]. The dynamics of Rayleigh-Duffing oscillator has been investigated 
widely in these years [Miwadinou,]. For example, in their work, Siewe Siewe 
et al. [2006] studied the nonlinear response and suppression of chaos by weak 
harmonic perturbation inside a triple well </> 6 —Rayleigh oscillator combined to 
parametric excitations. In [Siewe Siewe et al., 2009], the authors focussed their 
analysis on the occurrence of chaos in a parametrically driven extended Rayleigh 
oscillator with three-well potential modeled by: 

x — )i{l — x 2 )x + (1 — 7 coso :t){x + ax 3 + (3x 5 ) = F cos cot. (1) 

On the other hand, the Rayleigh-Duffing oscillator is used to model the roll 
motion of ships. The nonlinear ship rolling response generally is given by [Frances- 
cutto and Contento, 1999; Holappa et al., 1999; Scolan, 1999; Taylan, 2000] 

x + 2 /xx + (3\\x\x + / 32 x\x\ + S\x 2 x + 82 X 3 + +... + co^x + a?,x 3 + a$x b + 

... = Fcosixt, (2) 

where lj 0 and c 0 are internal and external frequencies respectively; //, and 8 t are 
linear, quadratic nonlinear and cubic nonlinear damping coefficients respectively 
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and Oj are restoring coefficients. F is the external excitation amplitude. In these 
cases, Francescutto and Contento [1999] used experimental results and parameter 
identification technique to study bifurcations in ship rolling modeled as : 
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x + 2fix + S 2 x 3 + UgX + cijxi = F cos c ot. (3) 

3 = 1 

Scolan [1999] applied the Melnikov method to nonlinear ship rolling in waves 
modeled by 

x + 2 fix — x(x 2 — x 2 )(x 2 — x 2 ) = Acosiojt + cfi). (4) 

The author provided that the heteroclinic orbits still exist whatever the “small¬ 
ness” of the perturbation as soon as the system is undamped. The existence of 
such cancellation is otherwise confirmed from an analysis of the erosion of the 
attraction basin. Application of the extended Melnikov’s method is used by Wan 
and Leigh [2008] for simple single-degree-of-freedom vessel roll motion. In their 
paper, the authors considered the two following ship rolling models: 

x + Six + e5 2 x\x\ + x — ax 3 = ef cos cut, (5) 


x + (3x + x(l — x)(l + ax) = F sin cot. ( 6 ) 

The authors compared their results with those obtained by Falzarano [1990] and 
Spyrou et al. [2002]. Our aim is to make a contribution to the study of the transi¬ 
tion to chaos in the modified Rayleigh-Duffing oscillator with three well potential 
possessing both homoclinic and heteroclinic orbit by using Melnikov’s theory, and 
then see how the fractal basin boundaries arise and are modified as the damping 
coefficient is varied. We also focus our attention on the numerical investigation of 
the strange attractor at parameter values close to the analytically predicted bifur¬ 
cation curves. 

The paper is organized as follows. In Section 2, we describe and analyze 
the model. Section 3 deals with the conditions of existence of Melnikov’s chaos 
resulting from the homoclinic and heteroclinic bifurcation. In Section 4, the basins 
of attraction of the initial conditions are droped to verify the effectiveness of the 
method. We provided a conclusion and summary in the last section. 

2 Description of the model 

We examine the dynamical transitions in periodically forced self-oscillating sys¬ 
tems containing the cubic and quintic terms in the restoring force; pure, hybrid 
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quadratic and cubic terms in nonlinear damping function as follows: 

x + n{l — x 2 )x + /3x 2 + k\xx + x + \x 3 + 5x 5 = F 0 cos t at, (7) 

where //, f3, hi, A, 8 , F 0 and c o are parameters. Physically, /j, /? and fci represent 
respectively pure cubic and pure, unpure quadratic nonlinear damping coefficient 
terms; oj and F 0 are respectively the frequency and the amplitude of external peri¬ 
odic forcing. Moreover A and 8 characterizes the intensity of the nonlinearity and 
e is the nonlinear damping control parameter. The nonlinear damping term cor¬ 
responds to the modified Rayleigh oscillator, while the nonlinear restoring force 
corresponds to the Duffing oscillator, hence the oscillateur is the so-called three- 
well potential modified Rayleigh-Duffing oscillator. 

We derive the fixed points and the phase portrait corresponding to the system 
(Eq. 0) when this is unperturbed. Eq. 0 becomes an unperturbed system and 
can be rewritten as 


x — y, y = —x — Xx 3 — 8x 5 , ( 8 ) 

which corresponds to an integrable Hamiltonian system with the potential function 
given by 


V(x) = \x 2 + \xx A + TT&c 6 , (9) 

2 4 b 

whose associated Hamiltonian function is 

H{x,y) = \y 2 + \x 2 F^Xx* F ^-Sx* 3 . (10) 

Z Z 4 D 

From Eqs. 0 and ( flQ| ), we can compute the fixed points and analyze their stabil¬ 
ities. 

• If A 2 < 4 8, the unperturbed system has only fixed point (0,0) which is a 
center and the potential V has single well (see Fig. |T|). 

• For A 2 = 4 8 with A > 0 and 8 < 0, the unperturbed system has three fixed 
points: two saddles connected by two heteroclinic orbits and one center. 

• For A 2 > 4 8 with A < 0 and 8 > 0, the unperturbed system has five fixed 
points: two saddles connected by two heteroclinic orbits, and the two saddle points 
are connected to themselves by one homoclinic orbit. The system has three centers 
and the potential defined by Eq. 0 has three-well (see Fig. [2]). 

In this paper, we consider the first and third cases. We set A = 0.3, 8 = —0.08 
and A = — 0 . 6 , 8 = 0.05 respectively in single-well potential and three-well 
potential cases. 
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Figure 1: (a) : Single-well potential of the unperturbed system; (6) : Shape of 
single-well potential for different values of S. 

3 Melnikov analysis 

In this section, we discuss the chaotic behavior of the system 

x + fi{ 1 — x 2 )x + /3x 2 + k\xx + x + Xx 3 + Sx 5 = F 0 cos u>t, (11) 

where y, hi,/3, /c 2 , A, u> and F 0 are assumed to be small. Hence, we write the 
dynamical system as 

x — y, y = —x — Xx 3 — Sx 5 — y(l — x 2 )x — /3x 2 H— k\xx + F 0 cos cot. ( 12 ) 

Classical Melnikov method is used in many cases to predict the occurrence of 
chaotic orbits in nonautonomous smooth one-degree-of-freedom nonlinear sys¬ 
tems [Litak, 2008; Sanjuan, 1999; Siewe Siewe, 2004; Siewe Siewe, 2005; Soli- 
man, 1992], It involves transverse intersection of stable and unstable manifolds 
that represent the starting point for a successive route to chaotic dynamics. Al¬ 
though a Melnikov theory is merely approximative, it is one of a few methods 
allowing analytical prediction of chaos occurrence. It enables prediction of values 
of the parameters associated only with the so-called heteroclinic or homoclinic 
chaos. This implies the existence of fractal basin boundaries, and so-called horse¬ 
shoes structure of chaos. To deal with such a question, we first derive the equation 
for the separatrix. 
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Figure 2: (a) : Triple well potential of the unperturbed system; ( b ) : Shape of 
triple well potential for different values of S. 


When the pertubations are added, the homoclinic or heteroclinic orbits might 
be broken transversely. And then, by the Smale-Birkoff Theorem [ Litak, 2008; 
Sanjuan, 1999; Scolan, 1999; Siewe Siewe et al., 2006], horseshose type chaotic 
dynamics may appear. It is well known, that the predictions for the appearance of 
chaos are limited and only valid for orbits starting at points sufficiently close to 
the separatrix. On the other hand it constitutes a first order perturbation method. 
Although chaos does not manifest itself in the form of permanent chaos, and some 
sorts of transient chaos may showup. However, it manifests itself in terms of 
fractal basin boundaries as it was shown in [Siewe Siewe et al., 2006], Let us 
start analysis from the unperturbed Hamiltonian Eq. ( flO] ). At the saddle point 
x = 0 , for an unperturbed system, the system velocity reaches zero so that the 
total energy has only its potential part. 

We apply the Melnikov method to the system in order to find the necessary 
criteria for the existence of homoclinic bifurcations and chaos. The Melnikov 
integral is defined as 

/ +oo 

/ ( x h, y h ) A g(xh, y h )dt, (13) 

-oo 

where the corresponding differential form / means the gradient of unperturbed 
hamiltonian while g is a perturbation from Eq. (fT2]). Eq. ([13]) can be rewritten as 
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follow: 


M ± (t 0 ) = ~l-i j y 2 h dt + p J yldt-ki J x h y 2 h dt - f3 J y* h dt+ 

~h J x h y\dt + F J y h cosuj(t + t 0 )dt, (14) 

where t 0 is the cross-section time of the Poincare map and t 0 can be interpreted as 
the initial time of the forcing term. 

Transforming Eqs. (|9|) and (10), for a closen nodal energy (H = 0) we get the 
following expression for velocity: 


dx I , 1 . A , S .. 
y = — = \ 2(— x 2 - x 4 - x 6 ). 

y dt V V 2 4 6 ’ 


Now one can perform the integration over x: 


t— to — ± 


dx 


\j~ l ~ ~ i x4 


(15) 


(16) 


X 


3.1 A single well potential case 

In this case, the fixed points are connected by an the heteroclinic trajectories given 
by [Siewe Siewe et al., 2004] 


Xh = ± 


cry sinh fit 


yh = ± 


a/ 2[1 + (1 - 7 2 ) coshf 2 (t - t 0 )] ’ 

a.'yfl cosh fit 

\/ 2 [l + (1 - 7 2 ) cosh fl{t - t 0 )] 3 ’ 


(17) 


where 


fl = x 2 p 


-(5(1 + p 2 


a 2 = xi(p 2 + 3y, 7 = 


2 P 2 


2 A + VA 2 - 4(5 
P = 


3(1 + P 2 )’ 


-A +VA 2 -4(5' 


Note that the central saddle point x 0 = 0 is reached in time t corresponding to 
Too and — oo respectively. 

After some algebra, the Melnikov function can thus be evaluated and a neces¬ 
sary condition for the onset of Melnikov chaos is given by 


F > 


a 3 fl 3 -\/ 2 (l — 7 2 ) sinh(|^) 
647to;7 


//(To I T,) 


8 « 3 7 3 r 2 2 / 5(7 — 67 s ) 
105(1 - 7 2 ) 3 / 2 


, (18) 
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where 


T 0 = 


a 2 n 

~T~ 


( 37 2 -l 
Vl-7 2 


+ 


3 7 2 + l 
2 7 


In 


l + 7\ 

1-7/ 


and 


Ti = 


315 7 8 —4207 6 +387 4 +207 2 +15 

( 1 — 7 2 ) 2 


g 4 Q 3 v 
I28O7 2 A 

15(3 7 2 +l)(l+7 7 4 ) ln 1+7 \ 
27 I-7 J ' 


Fi§ s UK 0 )’ (6) show the critical external forcing amplitude for different values 
of /1 and (3 respectively. The region under the curve represents the domain leading 
to suppression of horseshoes chaos. It appears that, depending on the parameters 
of the systems, the structure could or could not present chaotic dynamics. 



Figure 3: The critical curves for chaotic motions of a modified Rayleigh-Duffing 
oscillator in a single well potential case in (cu, F ) plane with A = 0.3, 6 = —0.08; 
(a) : effect of pure cubic damping coefficient /i and ( b ) : effect of pure quadratic 
damping coefficient (3. 


3.2 A three wells potential case 

In this subsection, the three wells type (the bounded double hump potential: Figj2]), 
the system Eq. ([8]) also has two hyperbolic fixed points and each point process two 
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different types of orbit: a heteroclinic orbit connecting the two saddle points an 
defined as [Siewe Siewe et al., 2005]: 


6x\\/2{l - jQsinh(f) \ 

+ coshOt) 3 J 

and a symmetric pair of homoclinic trajectories connected each point to itself 
given by 


(%hei Vhe) I i 


XiV / 2 cosh(^) 
>/(\k + cosh Qt) ’ 


(pChoi Vho) 


where 



a:i\/2sinh(y) 
y/(— 4) + cosh#t) ’ 


± 


Ox i v / 2(l _ v h) cosh(7^) 
2 \J (—'ll + cosh#£) 3 


,( 20 ) 


0 = - 1)2; a 2 = x 2 2 (p 2 + 3); 'k = §^; 

_ A+VA 2 -4 5 
0 — A+VA 2 —4(5 ’ 

and 


*2 = 

Let us first consider the case of heteroclinic orbit. After substituting the equa¬ 
tions of the heteroclinic orbits Xh and ?/;, given in Eq. ( [T9] ) into Eq. ( |T4] ), we 
calculate the Melnikov function and a sufficient condition for the appearance of 
chaos in sense of Smale is given by 


2(1 + 2xl0 2 (13-^)P' 

x\0 X he + 105(1-tt) 

where 


0 smh(^) 

27ro;a:i 


Xhe 


2 + \k + 


160(1+ ^) 3 


arcsin 'k + | 
\/l - 4> 2 


1 + 2 ^ + 


Oh^Ahe 

160(1 + \k) 3 


with 

„ _ 402^ 4 + 42314) 3 + 5063^ 2 + 4448^ + 976 

“ /ie “ 3('k - l) 


( 21 ) 
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and 


— 


20824' 4 + 3696^ 3 + 5889^ 2 + 405^ + 135 
3(4/ - 1) 


From each of these relations, the threshold values of the parameters for which 
the stable and unstable manifolds to the heteroclinic point intersect are obtained 
and the complicated behavior occurs whenever the above condition is realized. 
We study the chaotic threshold as a function of only the frequency parameter u. 
A typical plot is shown in Fig. |4} in which the critical heteroclinic bifurcation 
curves are plotted versus the frequency parameter ui. Figsg] (a), (b) show the 
critical external forcing amplitude for different values of /.( and j3 respectively. 
One can see (Fig|4] (a)) that when the value of the cubic damping parameter // 
increases, the thresholds of the critical values for heteroclinic bifurcation of the 
harmonic excitation F increase. The contrary effect is observed with the pure 
quadratic damping coefficient /3 (Figs{4](6)). We concluded that the parameters F 
and [i have the same effect on the critical value for chaotic motions which is the 
opposit one of parameter (3. 



Figure 4: The critical curves for heteroclinic chaos of a modified Rayleigh- 
Duffing oscillator in a single well potential case in (u. F ) plane with A = —0.6, 
S = 0.05; (a) : effect of pure cubic damping coefficient /i and ( b ) : effect of pure 
quadratic damping coefficient (3. 

Now, consider the case of homoclinic orbit. Substituting Eq. ( |20| ) into Eq. 
we calculate the Melnikov function. It is known that the intersections of the 
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homoclinic orbits are the necessary conditions for the existence of chaos. The 
Melnikov function theory measures the distance between the perturbed stable and 
unstable manifolds in the Poincare section. If M ± (fO) has a simple zero, then a 
homoclinic bifurcation occurs, signifying the possibility of chaotic behavior. This 
means that only necessary conditions for the appearance of strange attractors are 
obtained from the Poincare-Melnikov-Arnold analysis, and therefore one always 
has the chance of finding the sufficient conditions for the elimination of even 
transient chaos. Then the necessary condition for which the invariant manifolds 
intersect is given by 


F > 


e 


2x\ sin(#) 


2\[2k\ 


+ 


4(1+ $) 2 15(1 +^) 2 


( 22 ) 


where 


Xho = + 3^ - 2 - 


9 2 x\'E 


ho 


arcsin + — | 


120(1 + 4/) 2 

4t(34t - 1) - 02x2lAho 


x 


120 ( 1 + ^) 2 


with 


1 —ho 


54t 5 + 424' 4 + 6034' 3 + 561^ 2 + 566^ + 112 
4> - 1 ’ 


and 


A-ho = 


320^ 4 + 4004' 3 + 7654' 2 + 315^ + 90 
4/ - 1 ‘ 


These conditions provided a domain of the parameters where the system has 
transverse homoclinic orbits resulting in possible chaotic behavior. In order to 
have a visual information, we have plotted in Figs|5](a), (6) the dependence of the 
external excitation F on the frequency ui for different values of /i and k\ respec¬ 
tively. These figures, show that the parameters /a and k\ have the opposit effect on 
the homoclinic critical value for chaotic motions. 


4 Numerical simulations and analysis 

In order to verify the analytical results obtained in the previous sections, we have 
numerically integrated the system by using a fourth order Runge-Kutta algorithm 
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Figure 5: The critical curves for homoclinic chaos of a modified Rayleigh- 
Duffing oscillator in a triple well potential case in (c o, F ) plane with A = —0.6, 
6 = 0.05; (a) : effect of pure cubic damping coefficient [i and ( b ) : effect of 
unpure quadratic damping coefficient k\ 


to investigate the heteroclinic and the homoclinic chaos in the model. We con¬ 
sider the single well and triple well cases. We checked the effect of the nonlinear 
damping terms on the motion of the oscillator. In the same vein, we investigated 
how the basins of attraction are affected as the parameter /i, /?, k\ and F varied. 
We see through Figsj6}|7]and[8]that the basin boundaries become fractal or basin 
of attraction is destroyed. The blue zone stands for the area where the choice of 
the initial conditions leads to a chaotic motion while the white area is the domain 
of periodic or quasi-periodic oscillations. This means that the damping parameter 
H,/3,ki and amplitude of external forced F have contributed to These numerical 
results obtained confirm exactly the analytical results. 


5 Discussion and conclusions 

In this paper, the dynamics of a forced nonlinear damped system, the 0 6 modified 
Rayleigh-Duffing oscillator has been studied. The analytical criteria for the ap¬ 
pearance of chaos in the sense of Smale have also been derived using Melnikov 
theory. The effects of nonlinear damping and amplitude of external forced on 
Melnikov critical values have been investigated. In the single well potential case, 
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Figure 6: Basin of attraction in heteroclinic chaos of a modified Rayleigh-Duffing 
oscillator in a triple well potential: ( a)F = 0.05, (b)F = 0.1, (c)F = 0.4 and 
(d)F = 0.85, the others parameters are A = —0.6,5 = 0.05, k\ = 0.05,/? = 
0.05, /i = 0.0005, oj = 1. 

we noticed that the nonlinear damping parameters /i and /? and the forced ampli¬ 
tude F have the same effects. Melnikov critical value increases with each of these 
parameters. The parameters /.( and F have the same effects on homoclinic and het¬ 
eroclinic Melnikov criterion in three well potential case but the quadratic damp¬ 
ing f3 and k\ have the contrary effects. It should be noted that the pure quadratic 
damping affected the heteroclinic bifurcation and unpure quadratic damping af¬ 
fected the homoclinic bifurcation. A convenient demonstration of the use and 
accuracy of the method, is obtained from the basin of attraction. By means of 
the basin of attraction, we have shown that for certain regions of parameter space, 
the deterministic system driven harmonically experience behaviors that may be 
chaotic or non-chaotic. 

Finally, the above assessment reveals that roll damping is a very critical param¬ 
eter in motion characteristics of a ship. Therefore estimation of roll damping may 
lead to accurate prediction of chaotic ship rolling motions. In principle, it may be 
explained as the influence between nonlinear damping and excitation force. 
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Figure 7: Basin of attraction in homoclinic chaos of a modified Rayleigh-Duffing 
oscillator in a triple well potential: ( a)F = 0.1, (b)F = 0.4, (c)F = 0.5 and 
(d)F = 0.85, the others parameters are A = —0.6,5 = 0.05, k\ = 0 . 1 , f3 = 
0.05, n = 0.0005, u = 1. 
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Figure 8: Effects of nonlinear damping parameters on basin of attraction in the 
triple well potential (o) and ( b ) effect of F, (c) and (d) effects of k\ and (e) and 
(/) effects of ii. 





























